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Using the next-to-leading order low energy effective Hamiltonian for AB = 1 transitions, 
the effects of electroweak penguin operators in some two-body decay modes of B® meson 
are estimated in the Standard Model (SM). We find that in B Q S -> n + K~ and B° s -> K + K~ 
decay modes, the electroweak penguin effects are small, while in B° s — > n°K°, <f)K°, <fxp, 
the electroweak penguin operators enhance or reduce the pure QCD penguin and tree level 
contributions by 20% ~ 80% in decay width. We also present the results of CP asymmetries 
in these B° s deacy modes. 
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Recently Buras et.al. generalized the low energy effective Hamiltonian for \AB\ = 1 
transitions to next-to-leading order QCD corrections [1~4]. Now, we are in a position to 
investigate not only the pure QCD penguin and tree level contributions to B meson decays, 
but also the effects caused by electroweak penguin diagrams beyond the leading logarith- 
mic approximation. Naively, people believe that electroweak penguin contributions to B 
meson decays are suppressed by a factor of a em /a s ~ O(10~ 2 ) relative to QCD penguin 
contributions [5,6]. So electroweak corrections comparing with QCD penguin contributions 
in B system may be negligible. However this is not always true. In some B meson decay 
modes, electroweak penguin contributions can play a significant role [7,8]. In this paper, we 
investigate some B° s two-body decay channels. They are B° s — > -n + K~ , K + K~, ir°K , 00, 
and (j)K°, which are induced by QCD and electroweak penguin diagrams. We calculate the 
partial deacy widths, branching ratios, and the CP asymmetries of these decay modes. We 
compare the decay widths including full electroweak penguin contributions with that in- 
cluding only tree and QCD penguin diagrams. We find that in B® — > n + K" , and K + K~ , 
electroweak correction is small, it just modify the pure QCD corrections by 1% ~ 8%. While 
in Bg — > 7i°K°, (pK°, and 00, the electroweak penguin contributions are quite large relative 
to QCD penguin corrections. They enhance or reduce the results obtained by only taking 
into account tree and QCD penguin operators by 20% ~ 80%. The decay mode B® — > 00 
has been calculated in [8]. We include it here just for completeness and comparison. Our 
results for this channel agree with that in [8]. 

In the standard model, the number of colors N c is 3, while some people argue that the 
number favored by experimental data is 2 [9]. So in this paper we take not only N c = 3, but 
also N c = 2. We also present the results with N c = oo as in [10]. 

The next-to-leading order low energy effective Hamiltonian relevant to charmless B decays 
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can be taken as the following form [7] : 



H eff (AB = 1) = §6 



V2 



,q=u,c 



10 

k=3 



(1) 



where C fe (/x) (k=l,- • -,10) are Wilson coefficients which are calculated in renormalization 
group improved perturbation theory and include leading and next-to-leading order QCD 
corrections . v q is the product of Cabibo-Kobayashi-Maskawa (CKM) matrix elements and 
defined as 

VqdVqb for 6 ^ J transitions 



Vq = < 



V qs V* b for b — > s transitions 



In our numerical calculations the CKM matrix elements are taken as A = 0.22, A = 0.8, 
7] = 0.34, p = —0.12 in Wolfenstein parametrization, which are the preferred values obtained 
from the fit to the experimental data[ll]. The ten operators are taken as the following form 
[2,3] 

Q\ = {b a up) v ^ A {upq a )v-A, Q2 = (bu) V -A(uq)v-A, 

Qs = (bq) V - A J2(q'q')v-A, Qa = (bpq a ) V - A J2(q'aqp)v-A, 

q> q> 

Qb = (bq)v-A^2(q'q')v+A, Qe = (ppqa)v-A ^(q'^v+A, (2) 

q> q' 

3 - 

Qr =^W^E e ?'W)m 1 Qs = ^(bpq^v^AY.^'il'^v+A, 

<?' q' 

3 

Q9 = l(bq) v _ A J2 e q'(q'q')v-A, Qw = -{b^q a ) V -AY, e i^l' a qp)v-A^ 

q> q> 

where Q\ and Q\ are current-current operators, q=u, c, for q=c case, the two operators Q\ 
and Q2 are obtained through making substitution u — > c in Q\ and Q%. Q 3 ~ Q 6 are QCD 
penguin operators, the sum V] is runing over all the quark flavors being active at /x = m& 

q' 

scale, q'={u, d, s, c, b}. Qj ~ Qio are electroweak penguin operators, e q > are the electric 
charges of the relevant quarks in unit e which is the charge of the proton. The subscripts 
a, (3 are SU(3) C color indices. (V ± A) referes to 7^(1 ± 75). 



The Welsion coefficient functions Cj(/i) are renormalization scheme (RS) dependent be- 
yond the leading order approximation. Thus we should cancel this renormalization scheme 
dependence. So define renormalization scheme independent Wilsion Coefficients[12] 



C(/i) = (i + rja s (n)/4ir + fj a em (ii) / 4tt) ■ C(/j), 



(3) 



and treat the matrix elements to one- loop level [7], 



< Q T (/u) >=< Q T > 



(4) 



Combine eq.(3) and (4) we obtain 

< Q r (/i) • C(/x) > = <Q r >o-[i + ^K(^)-fJ T + ^(m e (/i)-f e ) T ]-C( /U ) 
= <Q T > -C'(^) 

(5) 

The corresponding elements of the matrices f s , r e , rh s (fj) and rh e ([/,) are given by [1, 4, 
13]. Substituting these matrix into eq.(5), we can obtain C[{ji) as 



a 



where P Sje are given by 



C' r = 



a = 



p. 



C 3 - P s /N c , 
C 5 - P s /N c , 
C 7 + P e , 

Cg + P e , 



C 2 = c 2 , 

C4 — C4 + Ps, 
Cq = Cq + P s , 

C's = C 8 , 
C' w = C10, 



(6) 
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G(m q ,q,n) 



3tt 



iVc 



y - G(m q ,q,ii) 



G(m g ,q,n) 



— 4 J x(l — x)dxln — - 



here q=u,c, for numerical calculation, we take m u = O.OOBGeV, m c = 1.35GeV, and use 
q 2 = ml/2, which represents the average "physical" value. 



We take m t = 174GeV,m b = B.OGeV, a s (M z ) = 0.118, a em (M z ) = and take the 
numerical values of the renormalization scheme independent Wilson coefficients Ci(/i) as [8] 

d = -0.313, C 2 = 1.150, C 3 = 0.017, 
C 4 = - 0.037, C 5 = 0.010, C 6 = - 0.046, 

(7) 

C 7 = - 0.001 -a em , C 8 = 0.049 ■ a em , C 9 = - 1.321 • a em , 
C 10 = 0.267 ■ a em , 
Using vacuum-saturation approximation the decay amplitude 

< XY\H e ff(AB = 1)|-B° > can be factorized into a product of two current matrix elements 

< X|J^|0 > and < F|J^|£?° >. The hadronic matrix elements are calculated in BSW 
method[14]. We define M££ w as in Ref.[10]. 

Ki<* = ^ < X\(q iq2 ) v . A \0 X Y\(bq 3 ) V - A \B >, (8) 
Now we can express the amplitudes < XY\H e ff(AB = 1)\B® > in terms of C[ and 



5 



M XY . 

<Jl<?2<?3 



<n+K-\H eff \B° s > 



<K+K-\H eff \B° s > 



< ir K°\H eff \B° s > 



< <M>\H eff \B° a > 



< <t>K*\H eff \B* > 



q=u,c 

Ml 



(m u + m d ){m b - m u 



v u a>i + v q {°3 



■(2a 5 + 3e M a 7 ) + a 9 • -e u } } ■ M^ U K , 



q=u,c 



M 2 3 
7- —Z J m -TT( 2a 5 + 3e M a 7 ) + a 9 • -e u } ] ■ M, 



K+K- 

USU 1 



+ 



_q=u,c 



q=u,c 
1 



3 3 



Mi 



°K° 



K + 1)( C i + CS-a,-a a .-e i+md{mi _ md) 



(a 5 + -e d • a 7 ) 



+(i + i)-M^ + cio)} 



M. 



ddd i 



q=u,c ^ iV c z 

+(l + i)(^ + Cio)-e s }]M/ s t 

Y V A a 3 + 



q=u,c 



■ - ■ r(-2a 5 - 3e s a 7 ) + -e s a 9 

m s + m d )(m s + m b ) 2 



+ 



Y «4 + «6 + ^e s (a 8 + aio) 

_q=u,c ^ ^ 



M, 



K°<j> 



sds 



(9) 



where the masses of d quark and s quark are taken as m d = O.OlGeV, m s = 0. 175(76 V" in 
numerical calculation, and is defined as 



a 2i-l 
d2i 



C" 



(i = 1,2, 3, 4, 5) 

The authers of Ref. [8] have calculated < (fxf)\H e ff\B® >. Our results in eq.(9) are con- 
sistent with theirs after making a charge conjugation. We follow the definition of relevant 
decay constants and form factors in Ref. [14], and use eq.(9~ 13) in Ref. [10] to calculate 
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\M*^ 2qs \ 2 . In the rest frame of B®, the two body decay width is calculated by 



q=u,c 



T(B° S -> 17) = —I < XY\H eff \B° s > | 2 ||, (10) 

where \p\ = K m b (m x +m y ) ){m b (m x m y ) )P . g ^ e ma g n i^ uc i e Q f momentum of the particles 
X and Y. M B , M x and M Y are the mass of B®, X and Y, respectively. For calculating 
the branching ratio, the total decay width of B° s is taken as Y tot = 4.91 x 10- 13 Gey[15]. 
We give our numerical results in Table 1-5, where "QCD+EW" means the decay width or 
branching ratio with full QCD penguin and EW penguin corrections, "QCD" with only 
QCD penguin corrections. r Q CD y^~ r Q CD represents the enhancement percentage of EW 
penguin contributions to the decay width. R ( ^^~~~ )) -^(f^)> anc ^ (^iTf ) represent 
the relevant ratioes of amplitudes of tree and penguins. 

From Table 1 and 2, we can see that for B® — > n + K~ and K + K~ , the EW penguin effects 
are small, they only reduce or enhance the decay width by 0.6% ~ 8%, and the results are 
not sensitive to N c . For B° s — > n + K~, the tree level contribution is dominant, while for 
B Q S — > K + K~, the QCD penguin contribution dominant. From Table 3 - 5, we find that 
for Bg — > n°K°, <pK° and 00, the EW penguin correction is significant. The percentage of 
correction is almost as high as —80.3% for B° s — > <f)K° decay modes when taking N c = 3. 

In Table 6 and 7, we present the branching ratios and CP asymmetry parameters of 
the five B° s decay modes discussed above. For / being non-CP-eigenstate cases, the CP 
asymmetry is defined as the following 

A T(B° a -> f) - r(jg -> f) nn 
cp r(a°-/) + r(5°-/)' [ } 

while for / being CP-eigenstate case, we calculate the CP asymmetry through[16] 

A = fo°°inB o phy s(t) -> /) - nB° phys (t) - mt 
cp jst^juo - /) + nK hys {t) - mt 

1 - |£| 2 - 2Im£(Am/T) 
~ (l + |e| 2 )[l + (Am/r)2]' 
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where Am/T ~ 20 which is the preferred value in the Standard Model [11]. 

We find that only for B° s — > 7r°i^°, and <fiK° two cases, CP violation parameter is sensitive 
to the number of colors N c , the other four decay modes are not. 

This work is supported in part by the China National Natural Science Foundation and 
the Grant of State Commision of Science and Technology of China. 
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Table 1 Numerical results of B° s — > ir + K 





T (10- 18 Ge\/) 


Fqcd+ew-Fqcd 

r QCD 


T) ( QCD+EW \ 
R ( tree ) 


R ( tree ) 




^UJJ+ri/ VV 




AT O 

l\ c — I 


1.83 


1.86 


-1 .Z.O /o 


n i Q7 
u. ry i 


U. lo ( 


U.U1U 


N c = 3 


2.016 


2.025 


-0.42% 


0.197 


0.193 


0.004 


N c = oo 


2.45 


2.42 


0.69% 


0.194 


0.200 


0.006 


Table 2 Numerical results of B° s -> 




T (10- 18 GeV) 


r QCD 


ry/ QCD+EW \ 

tree / 






QCD+EW 


QCD 


iV c = 2 


3.19 


2.93 


8.61% 


3.61 


3.43 


0.178 


iV c = 3 


3.58 


3.48 


2.93% 


3.62 


3.56 


0.065 


N c = oo 


4.22 


4.43 


-4.82% 


3.56 


3.68 


0.114 


Table 3 Numerical results of £° -> tt ^ 




T (10- 19 GeV) 


F QCD+EW qcd 


Y^t QCD+EW \ 
v tree / 


R( QCD } 

v tree ' 


v tree ' 


QCD+EW 


QCD 


iV c = 2 


1.59 


1.97 


-19.0% 


0.426 


0.617 


0.195 


iV c = 3 


0.45 


0.69 


-34.7% 


1.831 


2.516 


0.703 


iV c = oo 


0.59 


0.58 


1.72% 


0.501 


0.647 


0.149 



Table 4 Numerical results of 5° -> 





T (lO-^GeF) 


Fqcd+ew-Fqcd 
Tqcd 




QCD+EW 


QCD 


iV c = 2 


1.78 


3.46 


-48.6% 


iV c = 3 


0.12 


0.62 


-80.3% 


iV c = oo 


3.50 


2.22 


58.0% 



11 



Table 5 Numerical results of B° s — > 00 





T (10- 18 Ge\/) 


Fqcd+ew-Tqcd 

r QCD 




v^OLy-rii/ vv 


nrn 


N — 9 


2.46 


3.27 


9zL 8% 


iV c = 3 


1.56 


2.11 


-26.1% 


iV c = oo 


0.38 


0.52 


-26.9% 



Table 6 Branching ratioes and CP asymmetry parameters 
with only tree and QCD penguin contributions 



decay mode 


Br 


A 


iV c = 2 


iV c = 3 


iV c = oo 


N c = 2 


iV c = 3 


N c = oo 


B° s -> 7i + K- 


3.78 x 10~ 6 


4.12 x 10~ 6 


4.95 x 10~ 6 


-7.0% 


-7.9% 


-8.2% 


B° s -> X+K- 


5.97 x 10~ 6 


7.08 x 10~ 6 


9.02 x 10~ 6 


-2.31% 


-2.24% 


-2.18% 




4.00 x 10~ 7 


1.40 x 10~ 7 


1.18 x 10~ 7 


8.7% 


6.3% 


-32.6% 


S° -> (j)K° 


7.06 x 10~ 8 


1.25 x 10~ 8 


4.52 x 10~ 8 


-1.35% 


-5.0% 


0.97% 


B° s -> 00 


6.67 x 10- 6 


4.29 x 10~ 6 


1.06 x 10~ 6 


-0.010% 


-0.014% 


-0.024% 


Table 7 Branching ratioes and CP asymmetry parameters 
with full tree , QCD and EW penguin contributions 


decay mode 


Br 


A 

f^-cp 


N c = 2 


iV c = 3 


iV c = oo 


N c = 2 


iV c = 3 


N c = oo 




3.73 x 10~ 6 


4.11 x 10~ 6 


4.98 x 10~ 6 


-7.28% 


-8.01% 


-8.12% 




6.49 x 10~ 6 


7.29 x 10~ 6 


8.59 x 10- 6 


-2.22% 


-2.21% 


-2.24% 




3.24 x 10~ 7 


9.15 x 10~ 8 


1.20 x 10~ 7 


10.9% 


10.5% 


-32.0% 


- 0^o 


3.62 x 10~ 8 


2.48 x 10~ 9 


7.13 x 10" 8 


-1.58% 


-13.8% 


0.95% 


B° s - 00 


5.02 x 10~ 6 


3.19 x 10~ 6 


7.66 x 10~ 7 


-0.011% 


-0.016% 


-0.023% 
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